tsset 序号 年份

sum ln教育支出 ln医疗支出 ln消费品总额 ln土地增值税 坡度 贷款利率 湿地占比

gen ln教育支出=log(教育支出)
gen ln医疗支出=log(医疗支出)
gen ln消费品总额=log(消费品总额)
gen ln土地增值税=log(土地增值税)
gen ln收入=log(收入)
gen ln债务付息=log(债务付息)

reg ln销售价格  ln土地出让  ln赤字 lnGDP ln人口密度   三产比重 ln竣工面积  ln存款 

estat vif

reg ln销售价格  ln土地出让  ln赤字 lnGDP ln人口密度   三产比重 ln竣工面积  ln存款 
predict e, resid
gen abs_e=abs(e)
reg abs_e ln土地出让  ln赤字 lnGDP ln人口密度   三产比重 ln竣工面积  ln存款

sum ln销售价格 ln土地出让 ln赤字 lnGDP ln人口密度 三产比重 ln竣工面积  ln存款 ln土地增值税 贷款利率 ln债务付息 

pwcorr ln销售价格  ln土地出让  ln赤字 lnGDP ln人口密度   三产比重 ln竣工面积  ln存款, sig 

 
reg ln销售价格  ln土地出让,r  
reg ln销售价格  ln赤字, r  
reg ln销售价格  ln土地出让  ln赤字 lnGDP ln人口密度   三产比重 ln竣工面积  ln存款  i.序号  
reg ln销售价格  ln土地出让  ln赤字 lnGDP ln人口密度   三产比重 ln竣工面积  ln存款  i.序号   i.年份

//工具变量

reg ln销售价格  ln土地出让  ln赤字 lnGDP ln人口密度   三产比重 ln竣工面积  ln存款  i.序号   i.年份


ivregress 2sls ln销售价格  (ln土地出让=坡度 ln土地增值税)  ln赤字 lnGDP ln人口密度   三产比重 ln竣工面积  ln存款  i.序号, first


estat endogenous   //拒绝原假设，存在内生性
estat firststage  //弱工具变量检验
estat overid  //过度识别检验


gen 土利率=ln土地增值税*贷款利率
ivregress 2sls ln销售价格  (ln土地出让=土利率)  ln赤字  ln人口密度   三产比重 ln竣工面积  i.序号 i.年份, first

ivregress 2sls ln销售价格  (ln土地出让=土利率)  ln赤字 lnGDP ln人口密度   三产比重 ln竣工面积  ln存款  i.序号, first

estat endogenous   //拒绝原假设，存在内生性
estat firststage  //弱工具变量检验
estat overid  //过度识别检验



ivreg2 ln销售价格 ln赤字 lnGDP ln人口密度 三产比重 ln竣工面积  ln存款 (ln土地出让=土利率) i.序号


replace 债务付息=债务付息+1

ivregress 2sls ln销售价格  ln土地出让  (ln赤字=ln债务付息) lnGDP ln人口密度   三产比重 ln竣工面积  ln存款  i.序号 , first


estat endogenous   //拒绝原假设，存在内生性
estat firststage  //弱工具变量检验
estat overid  //过度识别检验

ivreg2 ln销售价格 ln土地出让 lnGDP ln人口密度 三产比重 ln竣工面积  ln存款 (ln赤字=ln债务付息) i.序号


ivregress 2sls ln销售价格   (ln土地出让 ln赤字 =土利率 ln债务付息) lnGDP ln人口密度   三产比重 ln竣工面积  ln存款  i.序号 , first

estat endogenous   //拒绝原假设，存在内生性
estat firststage  //弱工具变量检验
estat overid  //过度识别检验

ivreg2 ln销售价格   (ln土地出让 ln赤字 =土利率 ln债务付息) lnGDP ln人口密度  三产比重 ln竣工面积  ln存款  i.序号 



/////////////进一步分析///////////////
SFA2tier ln销售价格  ln土地出让  ln赤字 lnGDP ln人口密度   三产比重 ln竣工面积  ln存款 序号1-序号30, search   
   version 10.0
   ml max, difficult 

   SFA2tier_sigs 
   
 
   
   *---------------  
*-估计单边效应 
*---------------  
  drop   doc doc_pat
  SFA2tier_eff
  replace uw_diff_exp = -uw_diff_exp 
  
  foreach v of varlist u_hat w_hat uw_diff *_exp{
     replace `v' = `v'*100
  }
 
   *-样本总体
     tabstat u_hat_exp w_hat_exp uw_diff_exp,     ///
             s(mean sd p25 p50 p75) f(%6.2f) c(s) save
    tabstatmat all
    mat all = all'
    mat list all, title(full sample)



